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We demonstrate numerically and experimentally a conjugated gammadion chiral metamaterial that uniaxially 
exhibits huge optical activity and circular dichroism, and gives a negative refractive index. This chiral design 
provides smaller unit cell size and larger chirality compared with other published planar designs. Experi- 
ments are performed at GHz frequencies (around 6 GHz) and in good agreement with the numerical simulations. 



PACS numbers: 42.70.-a, 78.20.Ek, 42.70.Qs 



I. INTRODUCTION 

Materials whose magnetic/electric moment can be excited 
by the parallel external electric/magnetic field of the incident 
electromagnetic wave exhibiting optical activity^ are called 
chiral materials. They are characterized by the quantity of chi- 
rality, k = (tir — til)/2, where urIul is the refractive index 
of the right/left handed circular polarized wave (RCP/LCP). 
Natural chiral materials have very weak chirality (e.g., for 
quartz, ^~5xlO _5 atA = 400 nm). Five orders of magni- 
tude stronger chirality can be realized by chiral metamaterials 
made with sub- wavelength resonators. Chiral metamaterials 
recently attracted a lot of interest because of strong chirality, 2 
negative refractive index, 3,4 and the prospect of a repulsive 
Casimir forced Many chiral metamaterial designs have been 
proposed and demonstrated to obtain large optical activity, cir- 
cular dichroism, 6,7 and negative refractive index£~— 

In this article, we study both numerically and experimen- 
tally the conjugated gammadion chiral metamaterial that uni- 
axially exhibits huge optical activity, circular dichroism, and 
negative refractive index. This new design more clearly lacks 
any mirror symmetry plane compared to previous work*±2 In 
Ref. [13], the structure itself has a mirror symmetry plane. 
The lack of mirror symmetry only originates from the pres- 
ence of the substrate and different trace widths in the bi-layer 
structure due to fabrication constraints. While in our new 
design, the trace widths are the same and the lack of mir- 
ror symmetry is an inherent property of the particular design. 
And meanwhile, our new design possesses smaller unit cell 
size and larger chirality compared with other published planar 
designs: twisted-rosettes, 8 twisted-crosswiresr*^ and four-U- 

SRRs,ZJ±ti. 



II. EXPERIMENTAL AND SIMULATION 

The layout of the conjugated gammadion chiral metamate- 
rial is shown in Fig. 1. A 30 x 30 array of the conjugated 



gammadion resonator pairs is patterned on each side of the 
FR-4 board. The relative dielectric constant of the FR-4 board 
is e = 4.2 with a dielectric loss tangent of 0.02. The two lay- 
ers of the gammadion resonators are conjugatedly arranged in 
order to break the mirror symmetry along the direction per- 
pendicular to the metamaterial plane. The metamaterial pos- 
sesses C4 symmetry on the z-axis and, therefore, exhibits uni- 
axial chirality for the normal incident electromagnetic wave. 
The dimensions of the unit cell and the photograph of the ex- 
perimental sample are shown in Fig. Q] 




FIG. 1: (Color online) Scheme of the conjugated gammadion chiral 
metamaterial. Two gammadion copper resonators are conjugatedly 
aligned on each side of the FR-4 board. The geometry parameters 
are given as a x — a y — 10 mm, I = 8.1mm, w — 0.7 mm, 
d = 1.6 mm, and the thickness of copper is t = 0.036 mm. 



The amplitudes and phases of the linearly polarized trans- 
missions, T\\ (polarization of the transmitted wave is paral- 
lel with that of the incident wave) and Tj_ (polarization of 
the transmitted wave is perpendicular to that of the incident 
wave), and reflection coefficient R (R = Ru, R± = 0)^ 
are measured using an HP 8364B network analyzer with two 
Narda standard horn antennas. The circular polarized trans- 
missions of the RCP and LCP, Tr and T^, and reflections, Rl 
and Rr, can be obtained from the linearly polarized transmis- 
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sion coefficients by: 

T R = T,| - iT^T L = T,| + zT ± , Rr = Rl = R. (1) 

The circular polarized waves are the eigenwave functions of 
the chiral metamaterials. The cross transmissions, from LCP 
to RCP and vice versa, are zero. The numerical simula- 
tions were performed using the frequency domain solver of 
the CST Microwave Studio (Computer Simulation Technol- 
ogy GmbH, Darmstadt, Germany), which implements a finite 
element method. In the simulations, the unit cell boundary 
condition was applied and the circular polarized eigenwaves 
used directly. 



III. RESULTS 

Figures. [2a) and (b) show the simulation (left) and experi- 
mental (right) results of the transmissions and reflections. The 
calculated results nicely agree with our experimental results. 
There are two resonances in the transmission spectrum. One 
corresponds to a peak around / = 5.6 GHz and the other, a 
dip, around / = 7.8 GHz. The first resonance is much sharper 
than the second one. The feature of the transmission peak for 
both RCP and LCP at the first resonance is observed for the 
first time in chiral metamaterials. From the following retrieval 
results, we conclude that this peak is induced by the strong 
magnetic response, which decreases the degree of impedance 
mismatch. At the resonances, the transmission spectra for the 
RCP and LCP waves are significantly different. At the first 
resonance, the transmission of LCP is larger than that of the 
RCP; at the second resonance, it reverses. This difference be- 
tween the amplitudes of two transmissions is characterized by 

ellipticity, ry = - tan -1 ( |^| 2 ~|^| 2 ), shown in Figs. [2te) and 
(f). Far from the resonances, the phases of the RCP and LCP 
converge. In the vicinity of the resonances, they are obviously 
different, which will induce the rotation of the polarization 
plane of a linearly polarized light as it passes through the chi- 
ral metamaterial. The difference between the phases is char- 
acterized by the rotary power, = |[arg(T^) — arg(T#)], 
shown in Figs. Oc) and(d). The most interesting region is 
around r] = 0, which will give the "pure" optical activity ef- 
fect, i.e., for the linearly polarized incident wave, the trans- 
mission wave will still be linearly polarized just with the po- 
larization plane rotated by 6 degree. This design possesses 
very strong optical activity. At r] = (i.e., / = 6.55 GHz), 
the rotation angle is about 30 degrees for 1.6 mm thick meta- 
material, i.e., the rotation angle is as large as 856° /A. 

After obtaining the data of the transmissions and reflec- 
tions, we can apply a retrieval procedure to obtain the effective 
constitutive parameters. 17,18 The chirality k can be obtained 
directly from the transmissions as 

arg(T L ) - arg(T^) + 2m?r 
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FIG. 2: (Color online) The simulation (left) and experimental (right) 
results of the transmissions and reflections, and the rotation angle 
and circular dichroism. 



Simulation 



Experiment 



c 

—i 
c 



10 
5 


-5 

20 
10 

-10 

150 
100 
50 

-50 



. — — / 
/ 


— n 

K / 


N 




_ " ft 

(M. 




(a) 


' " " N 

/ 


1 
1 


'(b) 






/' — n L 

i 1 / 
/ 1 — n R / 

' * ' 1 








1 _ _ - " 






X "T* 

■/ 

(c) ' 




■(d) ' 




"(e) " 




1 
i 
/ 


-"- ^ 1 




(0 





5 6 7 8 
Frequency (GHz) 



5 6 7 8 
Frequency (GHz) 



FIG. 3: (Color online) The retrieval results from the numerical (left) 
and experimental (right) results of the transmissions and reflections. 



where fco is the wave vector in the vacuum; d is the thick- 
ness of the sample (here, the thickness of the retrieved sam- 
ple is chosen as 1.6 mm, i.e., the thickness of the FR-4 
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board) and the integer, m, is determined by the condition of 
— 7r < arg(Tk) — arg(Tft) + 2mir < tt for one unit cell. In 
previous studies the retrieved parameters depend on the 
thickness of the sample 19 and once we have multi-layer sam- 
ple, the retrieved parameters, for the weakly coupled system, 
are very close to the one unit cell and converge very fasA 
For the strongly coupled system, the retrieved parameters are 
completely different from the one unit cell system.^^I 

Comparing with Figs. Etc-f), we note that the real and 
imaginary parts of k relate to the azimuth rotation angle, 0, 
and the circular dichromism, 77, respectively. The average re- 
fractive index, n = (n^ + ur)/2, and the impedance z can 
be obtained by using the traditional retrieval procedure 18 after 
taking the geometric average transmission, T = ^/TlTr. The 
branch of the square root is chosen coordinating with the real 
part of k. Then, the other parameters can be calculated by: 
n± = n ± k, e = n/z, and /i = nz. Figure [3] shows the re- 
trieval results. The chirality is very large, n = 2.35 at r\ = 0, 
which corresponds to frequency, / = 6.55 GHz (k = 2.09 for 
the experimental result). This strong chirality can easily push 
the refractive index of the RCP/LCP eigenwave, n± = n ± k, 
to become negative at the resonances as shown in Fig. Oc) 
and (d). Due to the loss of the PC board, the figure of merit 
(FOM), — Re(n)/Im(n), is very low compared with the tra- 
ditional negative index metamaterials.— If using a lower loss 
substrate, the FOM can be larger than 10 according to our 
numerical simulation. In Figs. [3e) and (f), we show the re- 
trieved real parts of the permittivity and permeability of the 
chiral metamaterials. 



IV. DISCUSSION 

To illustrate the origin of the chiral response of our meta- 
material, in Fig. 01 we discuss a procedure of transmuta- 
tion from the simple ^-particle chiral element to the conju- 
gated gammadion chiral metamaterial. The ^-particle chiral 
element is one of the important chiral structures studied an- 
alytically elsewhere. 1 Here, we place two ^-particle chiral 
elements together to form a conjugated ^-particle chiral el- 
ement pair. The linearly polarized external electric field can 
drive the electric dipole via the vertical arms. This electric 
response (ER) generates circular electric current on the loop 
which gives a magnetic moment in the same direction as the 
electric dipole. Rotating the conjugated ^-particle chiral ele- 
ment pair by 90 degrees, the external magnetic field will drive 
the magnetic moment via the loop. This magnetic response 
(MR) then drives the current on the vertical arms which gives 
an electric dipole in the same direction as the magnetic mo- 
ment. This clearly shows how the magnetic/electric moment 
is induced by the electric/magnetic field (of the incident EM 
wave) in the parallel direction. According to the aforemen- 
tioned definition, this results in the chirality. If adding the 
two ^-particle pairs together, it will be a bi-isotropic uniaxial 
chiral metamaterial. After stretching the semi-circle loop to a 
straight line, it becomes to our current planar design. 

To further illustrate the origin of chirality in our conju- 
gated gammadion chiral metamaterials, Fig. \5\ shows three 




FIG. 4: (Color online) The procedure of the transmutation from the 
simple Q-particle chiral element to the conjugated gammadion chiral 
metamaterial. The procedure shows how the magnetic/electric mo- 
ment is induced by the electric/magnetic field (of the incident EM 
wave) in the parallel direction. ER and MR are the electric and mag- 
netic response modes, respectively. The arrow indicates the current 
direction. 




FIG 5: (Color online) Three basic electromagnetic responses in con- 
jugated gammadion chiral metamaterials. The arrow indicates the 
current direction. ER and MR are comparable with those in Fig. [4] 
ER' is a new emerged electric response which is driven by the central 
arms. ER' can not give chirality response. The incident wave is the 
same as in Fig. 0] 
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FIG. 6: (Color online) Illustration of superposition at / = 5.6 GHz 
(upper) and / = 7.8 GHz (down). 
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basic electromagnetic responses (ER, MR, and ER') in con- 
jugated gammadion chiral metamaterials. ER and MR are 
comparable with those in Fig. [4] ER 7 is an additional, non- 
resonant electric response which is driven by the central par- 
allel arms. ER' can not give chirality response because there 
is no magnetic moment generated. In our conjugated gamma- 
dion chiral metamaterials, each resonance is a superposition 
of these three electromagnetic response modes. For example, 
as shown in Fig. [6] the resonance at f=5.6 GHz can be consid- 
ered as the superposition of ER, MR, and ER 7 . The existence 
of the MR is also reflected by the retrieval results in Figs. Oe) 
and Of), where the effective permeability has a strong reso- 
nance at f=5.6 GHz. And the resonance at f=7.8 GHz can be 
considered as the superposition of ER and ER 7 , then the re- 
trieval results only show the electric response. 



ure[7]shows that the retrieved refractive indices for the multi- 
layer metamaterial is rapidly converged. Three layers (N=3) 
almost overlaps with two layers (N=2). Therefore, the sin- 
gle layer we studied previously can represent the properties 
of bulk chiral metamaterials for these weakly coupling cases. 
The differences between the refractive indices in Fig. [3] and 
Fig. [7] for the single unit cell are due to the different thick- 
nesses of the unit cell. In Fig. [3l we take the unit cell thick- 
ness as 1.6 mm, while it's 6.4 mm in Fig. [7J Comparing Fig. 
[7] and Figs. [3te) and[3f), the magnitude of the retrieved ef- 
fective parameters of ur and n l decrease, as the size of the 
unit cell increases^ Note, for the strong coupling cases (i.e., 
small D), the strong coupling between each layer will induce 
differences between one layer and multi-layer structure 21 and 
slow the convergence of the retrieved parameters i22 
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FIG. 7: (Color online) The retrieved refractive indices for multi-layer 
conjugated gammadion chiral metamaterial 



V. CONCLUSION 

In conclusion, we have designed and studied a conjugated 
gammadion chiral metamaterial at around 6 GHz. This chiral 
metamaterial exhibits huge uniaxial optical activity. The ro- 
tation angle is as large as 30° for 1.6 mm thick metamaterial, 
corresponding to k = 2.35 with r] = 0. The chiral metamate- 
rial also exhibits uniaxial negative refractive index, due to its 
strong chirality. The origin of the chirality is intuitionally and 
physically explained, based on the ^-particle chiral element 
model. 



Some may raise the question whether our design is just a 
thin meta- surface that exhibits rotation of polarization phase 
or can represent the properties of the bulk chiral metamaterial 
as we expected. In order to clarify this point, in Fig. Owe plot 
the retrieved refractive indices for multi-layer metamaterial (a 
pair of doule-layered conjugated gammadions is called one 
layer.) . Here, we study the weakly coupling case:— The sep- 
aration D between each layer of conjugated gammadion pairs 
is three times of the thickness of FR-4 board, i.e., D = 3d and 
the total thickness of each unit cell is D + d = 6.4mm. Fig- 
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